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Abstract 


Global climate changes and population growth have increased the need to exploit underground 
water resources. The phenomenon of land subsidence as a result of groundwater extraction and 
water depletion is one of the environmental hazards that threatens societies all over the world. In 
this research, for the first time, Sentinel-1 radar satellite images were used by the 
SNAP2StaMPS to automate the process of processing and analyzing the subsidence in Marand 
Plain. In recent years, the uneven development of agricultural lands and the excessive extraction 
of underground water in the Marand plain in East Azerbaijan province have caused landslides 
on the surface of this plain. In this study, an automatic interferometric processing algorithm 
between SNAP and StaMPS was used to determine the displacement rate in the direction of the 
satellite line of sight. The results of the PSI interferometry technique were performed on 133 
images of the Sentinel-1 satellite in the descending orbit of the Sentinel-1 satellite. Based on the 
obtained results, the annual land displacement rate for Marand plain in 2016, 2017, 2018, 2019, 
and 2020 was -13.7, -12, -2.15, -12.3 and -13.1 cm respectively. In order to validate the results, 
the amount of subsidence resulting from the interferometric processing was compared with the 
amount of groundwater level drop in the study area, using the unit hydrograph analysis method. 
The results showed that the amount of water level drop and the results of hydrograph analysis of 
observation wells in the region are consistent with the results of time series maps obtained from 
interferometry. Therefore, the automatic method presented in this study can be used for 
monitoring land subsidence. 
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جغرافیا و مخاطرات محیطی, سال یازدهم. شمار: چهل و چهارم زستان ۶۰۱ صص ۲۱-۶۲ 
مقاله پژوهشی 
ویژ‌نامه (چالش جهانی فرونشست زمین: مدیریت بحران یا بحران مدیریت) 
پایش پدیده فرونشست دشت‌ها بر مبنای الگوریتم خود کار SNAP2STAMPS‏ به روش تداخل سنجی 


راداری (PSD‏ (مطالعه موردی: دشت مرند) 
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حکیده 


در این تحقیق برای اولین بار به‌منظور خودکار سازی فرآیند پردازش و تحلیل فرونشست دشت مرند از 
تصاویر ماهواره‌ای راداری ستتینل-۱ توس ط SNAP2StaMPS c‏ استفاده‌شده است. در سال‌های 


اخیر توسعه ناهمگون اراضی کشاورزی و برداشت بی‌رویه از آب‌های زیرزمینی دشت مرند در استان 


آذربایجان شرقی موجب بروز پدید B39‏ = زمین در سطح این دشت شده است. لوله زایی میله 
چاه‌های منطقّه. شکاف‌های افقی در سطح این دشست. حکایت از رخداد فرونشست سطح زمین دارد. در 
این تحقیق برای تعیین نرخ جابجایی در جهت خط دید ماهواره. از الگوریتم پردازشی تداخل سسنجی 
خودکار بین نرم‌افزار StaMPS ;SNAP‏ استفاده شد. نتایج فن تداخل سنجی PSI‏ بروی ۱۳۳ a‏ 


سنتینل-۱ در مدار Descending‏ ماهواره 61001061-1 حاکی از نرخ جابجایی سالانه زمین برای دشت 


مرند از سال‌های ۳۰۱١‏ ۳۰۱۷ ۲۰۱۸ ۲۰۱۹ و ۲۰۲۰ به ترتیب eV‏ ۱۲ے Oe aA‏ - و ۱۳/۱ 
میزان افت سطح آب‌های زیرزمینی در محدوده مطالعاتی از طریق روش تحلیل هیدروگراف واحد انجام 


۱ نویسنده مسئول: ۰۹۳۸۷۲۰۷۹ Email: Samira.najafvand @ tabrizu.ac.ir‏ 
نحوه ارجاع به این مقاله: 
روستایی. شهرام؛ نجف‌وند. سمیرا؛ Ate)‏ پایش ody‏ فرونشست دشت‌ها بر مبنای الگوریتم خودکار SNAP2STAMPS‏ 


به روش تداخل سنجی راداری (PSI)‏ (مطالعه موردی: دشت مرند). جغرافیا و مخاطرات محیطی.4(۱۱),. صص ۲۱-۶۲ 
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Lal > YY‏ و مخاطرات محیطی شمارة چهارم 


شد. نتایج نشان می‌دهد که میزان افت سطح آب و نتایج حاصل از تحلیل هیدروگراف چاه‌های مشاهده‌ای 
منطقه با نتایج حاصل نقشه‌های سری زمانی حاصل از تداخل سنجی تطابق دارند؛ بنابراین» از روش 


خودکار ارائه‌شده در این پژوهش می‌توان جهت پایش فرونشست زمین استفاده کرد. 


کلیدواژه‌ها: "C‏ زمین روش تداحل oe‏ راداری «(PSD‏ دت U ya‏ الگوریتم خودکار 


Snap2stamps 


۱- مقدمه 


مخاطره فرونشست زمین شامل فروریزش یا نشست رو به پایین سطح زمین است. این پدیده به یکی از مسائل و 
مشکلات اساسی جوامع بشری به دلیل خسارات جبران‌ناپذیر نسبت به دیگر مخاطرات تبدیل‌شده است. GS JSS‏ 
این پدیده ب‌صورت آرام در طی مدت‌زمان طولانی اتفاق می‌افتد و فاقد برگشت‌پذیری است. فرونشست به‌طورکلی 
ناشی از دلایلی همچون انحلال تشکیلات زیرسطحی, تراکم رسوبات در اثر استخراج سیالات و ذخایر زیرزمینی» 
بارگذاری» زهکشی و یا ارتعاش تراکم هیدرولیکی. ناپایداری خاک‌ها در سطوح شیب‌دار انقباض خاک‌های رسی» 
ذوب بخ لایه‌های منجمد در اعماق زمین و ... پدیدار می‌شود. فرورفتگی قیفی شکل زمین» پس‌روی dos el‏ ریزش 
سطح و جابجایی پوسته‌ای زمین از علائم OF‏ است. علاوه بر این» فرونشست زمین در مناطق شهری می‌تواند باعث 
فروریختن فونداسیون تخریب خط لوله زیرزمینی و سیب جاده شود. این مخاطره موجب تغییرات توپوگرافیکی 
سطح زمین نیز می‌شود که در اکثر اوقات آسیب‌های جبران‌ناپذیری به بار می‌آورد. 

fle‏ وقوع فرونشسست زمین در نقاط مختلف دنا بهدلایل مختلفی رخ مدهت ازجملهء احیای مناطق ساحلی» 
توسعه فعالیت‌های صنعتی و ساخت‌وسازهای شهری در مناطق دلتایی و استخراج بی‌حد GET‏ زیرزمینی در 
جلگه‌ه؛ که جابجایی سطح زمین را به دنبال دارد. در سراسر جهان, دشت‌های ساحلی و مناطق دلتایی رودخانه‌ای 
بیشترین فرونشست زمین رادارند (هررا گارسیا" و همکاران,۲۰۲۱). 

محرکان اصلی فرونشست به دو عامل انسانی و طبیعی تقسیم می‌شوند» عوامل ناشی از فعالیت‌های انسانی شامل 
استخراج سیالات از لایه‌های زیرسطحی (مانند آب نفت و GIS‏ تغییر کاربری زمین (به‌عنوانمثال بارگذاری ناشی از 
ساختمان‌های ساخت بشر)» ساخت تأسیسات زیرزمینی و معدنکاری و فعالیت‌های طبیعی شامل گسل» حرکات 
زمین ماک کیت و فشردگی رسویانته فروریزش cond‏ افزایش سطح dios OT‏ اکسیداسیون و ژهکشی خاک‌های 
آلی» فرسایش کارستی و فروافتادگی ناشی از پرمافرست هستند (کروستو" و همکاران۲۰۱). در این میان» عامل 


] Herrera-García 
2 Crosetto 


سال یازدهم پایش پدیده فرونشست دشت‌ها بر مبنای الگوریتم خودکار .... ۳۳ 


زیرزمینی» به دلیل توسعه شهری و کشاورزی. محرک‌های اصلی فرونشست زمین. به‌ویژه در مناطق خحشک و 
نیمه‌عشک هستند (باقری گاو کش " و همکاران» ۲۰۲۱). 

امروزه فرونشست زمین یک مخاطره جهانی ناشی از عامل انسانی و طبیعی است و بیشترین مقدار میانگین 
فرونشست ساانه (LSavg)‏ در شرق چین» دشت مرند در شمال غرب ایران» ساحل خلیج تگزاس در جنوب 
ایالات‌متحده جاکارتای شمالی در اندونزی. منطقه شهری مکزیکوسیتی در مرکز مکزیک مشاهده‌شده است؛ بنابراین 
پیش‌بینی» تشخیص و کاهش فرونشست زمین در دشت‌ها از اهمیت بالایی برخوردار است (هررا گارسیا و همکاران, 
(vey)‏ 


در مطالعه حاضر جهت پایش فرو 


زمین در دشت مرند در استان آذربایجان شرقی از روش پراکنش 
گرهای دائمی PSD‏ که برای اولین بار توسط فرتی " و همکاران در سال ۲۰۰۱ مطرح شد استفاده‌شده است (فرتی و 
همکاران. ۲۰۰۱). این تکنیک مزایای زیادی دارد. مانند پوشش مکانی وسیع. حساسیت زياد به تغییر شکل‌های 
کوچک و عملکرد در روز و شب و تمام شرایط آب و هوایی را ارائه می‌دهد. که این مزایا پایش دقیق فرونشست 
زمین را تضمین می‌کند. دقت GUL‏ روش تداخحل سنجی راداری PSI‏ سبب شده است تا در بسیاری از تحقیقات از 
این روش به‌منظور برآورد میزان فرونشست استفاده شود. 

در ایران نیز بحث فرونشست زمین به یک بحران ملی تبدیل‌شده که در بیشتر موارد ناشی از برداشت بی‌رویه از 
منابع آب‌های زیرزمینی است. در این خصوص پژوهشگران مختلفی میزان فرونشست زمین را در مناطق مختلف 
ایران انجام‌شده است. 

خرمی YS‏ در فاصله زمانی سال ۲۰۱۶ تا ۲۰۱۷ میلادی با استفاده از 74 تصویر راداری مدارهای Ub‏ گذر و 
پایین گذر ماهواره ستتینل-۱) بیش‌ترین مقدار فرونشست را در شمال غرب مشهد و در حدود ۱۶ سانتی‌متر در سال 
به دست آورده است. مژمنی (۱۳۹۷) نرخ فرونشست را در دشت کبودآهنگ-فامنین همدان را در بازه‌های زمانی 
۳ ۲۰۱۰ و ۲۰۱۵ تا ۲۰۱۸ با روش پراکنش گرهای دائمی» به ترتیب -٩۲‏ و AA‏ میلی‌متر در سال محاسبه 
کرده است. فروغ نیا و همکاران (۱۳۹۷) با به‌کارگیری تصاویر سنجنده های انویست و ستتینل-۱ نرخ فرونشست را 
از سال ۲۰۰۳ تا ۲۰۱۹ حدود ۱۳۰ میلی‌متر در سال برآورد کرده‌اند. شرافت و همکاران (۱۳۹۸) در بازه زمانی ۲۰۰۳ 
تا ۲۰۰۵ با استفاده از ۱۲ تصاویر ماهواره‌ای Envisat‏ متوسط سرعت فرونشست در دشت oS pl‏ یزد را ۲ تا ۳ 


سانتی‌متر و در منطقه مهردشت ۵ سانتی‌متر به دست آورده‌اند. 


1 Bagheri-Gavkosh 
2 Ferretti 


yé‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


احمدی و همکاران (۱۳۹۹) نرخ فرونشست دشت اسداآباد همدان را ۲۹۱-۲۲۰ میلی‌متر در سال به دست 


آورده است. همچنین با توجه به نتایج حاصل از تحقیق نرخ فرونشست از مناطق شهری به سمت دشت‌های اطراف 
صعودی بوده است. پاپی و همکاران (۱۳۹۹) به تحلیل سری زمانی فرونشست در غرب تهران (دشت شهریار) و 


ارتباط oT‏ با برداشت آب‌های زیرزمینی با تکنیک تداحل سنجی راداری پرداخته‌اند. نتایج این تحقیق Sle‏ پیوستگی 


رخداد فرونشست در زمین‌های کشاورزی با نرخ متوسط Ae‏ سانتی‌متر و حداکثر ۲۳- سانتی‌متر در سال در منطقه 
دشتی نشان می‌دهد. که افت ۰/۵ تا ۱/۵ متری سطح آب در آبخوان علت آن است. خوش لهجه آذر (۱۳۹۹» میزان 
جابجایی بخشی از دشت‌های کبودرآهنگ- فامنین و همچنین رفتار زمین قبل از رخداد فرو چاله در کرداباد همدان 
را با دو روش PSI‏ و SBAS‏ به ترتیب ۵/۷۱ سانتی‌متر و ۲۶/۷ سانتی‌متر در سال به دست آورده است. قره چلو و 
همکاران (۱8۰۰)» با استفاده از تصاویر راداری آلوس (ALOS)‏ و ستتینل-۱ میزان فرونشست سالانه را در دشت 
مشهد بررسی کرده‌اند. نتایج حاصل از فرآیند تداخل‌سنجی راداری برای ماهواره آلوس در این تحقیق, بیانگر مقدار 
فرونشست ۲/۵ تا ۸/۳ و ۷/۶ سانتی متر برای تصاویر ALOS‏ و برای ماهواره سنتینل-۱ بین ۱/۱۳ تا ۶/۱۷ سانتی‌متر 
و ۳/۲۰ سانتی‌متر در بازه زمانی ۲۰۰۷ تا ۲۰۰۸ بوده است که نشان‌دهنده روند افزایش پدیده فرونشست در منطقه 
بوده است. 

ماتئوس ' و همکاران (۲۰۱۷)» فرونشست زمین در آبخوانی در اسپانیا را با استفاده از داده‌های راداری انویست در 
بازه زمانی ۲۰۰۳ تا ۲۰۰۹ و ماهواره (COSMO-SkyMed)‏ بین سال‌های ۲۰۱۱ تا ۲۰۱۶ و تصاویر سنجنده سنتینل- 
۱ بین سال‌های ۲۰۱۵ Y VG‏ مقدار فرونشست را 00 میلی‌متر به دست آورده‌اند. آن‌ها پیش‌ترین میزان فرونشست 
زمین مربوط به یک دوره طولانی و خشک در منطقه بوده که بر روی زمین‌های رسی رخ‌داده است. سان" و همکاران 
OOV)‏ طی پژوهشی فرونشست زمین در بخش‌های جنوبی پایین‌دست دشت لبائوهه " چین بین سال‌های ۲۰۰۷ تا 
۱ با استفاده از روش تداخل سنجی پراکنده ساز دائمی (PST)‏ پرداخته‌اند. تصاویر مورد پردازش در این تحقیق, 


تصاویر راداری مربوط به ماهواره L‏ باند آلوس پالسار بوده که از طریق نرم‌افزار استمپس StaMPS)‏ مورد پردازش 


قرارگرفته‌اند؛ و فرونشست سالانه زمین را در این دشت به دست آورده‌اند و نقشه‌های سری زمانی برای سال‌های 
۷ تا ۲۰۱۱ ایجاد کرده‌اند. در کل بر اساس نتایج روش پراکنده سازهای دائمی» سه ناحیه فرونشست در زیرزمین 
کشف شد و این مناطق در مناطق نفتی» مناطق استخراج نمک و مناطق ساحلی قرار داشتند. میزان فرونشست زمین 
از ۵۰ میلی‌متر در سال به ۲۳۱ میلی‌متر در سال رسیده است» مقصودی و همکاران YA)‏ نرخ فرونشست جاوا در 


غرب اندونزی را در دو 69b‏ زمانی ۲۰۰۷ تا ۲۰۰۹ و ۲۰۱۵ تا ۲۰۱۵ را طبق داده‌های ستتیئل-۱ و آلوس به ترتیب ٤/١‏ 


1 Rosa María Mateos 
2 Sun 
3 Liaohe Plain 


سال یازدهم پایش پدیده فرونشست دشت‌ها بر مبنای الگوریتم خودکار .... Yo‏ 


و ۲/۸ میلی‌متر در سال برای این منطقه گزارش کرده است. لیوسیس" و همکاران (۲۰۱۸) فرونشست در نواحی 
روستایی الویگان ‏ امارات را بین سال‌های ۲۰۱۰-۲۰۰۳ به ترتیب برای داده‌های آلوس و ستتینل-۱ حدود VA‏ 


سانتی‌متر و حدود Y‏ سانتی‌متر در سال به دست آوردند, دالیو "و همکاران us 23 9 (Y YA)‏ زمین را در 


باتلاق‌های نمکی لاگونی ونیز با استفاده تداخحل سنجی راداری و تحقیقات زمینی را موردبررسی قرار داده‌اند. این 
محققین از تکنیک پراکنده ساز دائمی (PSI)‏ داده‌های راداری در یک دوره ۵ ساله مربوط به منطقه موردمطالعه 
استفاده کرده‌اند. نرخ فرونشست را بیش از ۲۰ میلی‌متر در سال به دست آورده‌اند. عامل اصلی فرونشست عامل 
انسانی بوده است. میلر* و همکار (۲۰۱۹)» در مقاله‌ای به بررسی فرونشست زمین در هوستون و ارتباط آن با سیل 
حاصل از طوفان هاروی پرداخته‌اند. این محققین با استفاده از تکنیک تداخل سنجی راداری و تصاویر سنجنده 


سنتینل A=‏ سری B‏ نقشه پهنه‌های خطر سیلاب را مشخص کرده‌اند. نقشه فرونشست زمین را نیز قبل از وقوع طوفان 
را با استفاده از پردازش داده‌های چند زمانه SAR‏ ماهواره‌ای پیشرفته (ALOS)‏ و ماهواره‌های سنتینل-۱ سری B‏ به 
دست آورده‌اند. آن‌ها دریافتند که مقدار فرونشست GUL‏ £4 میلی‌متر در سال و ۳۶ میلی‌متر در سال در طی بازه زمانی 
پردازش تصاویر آلوس (ژوئیه ۲۰۰۷-ژانویه ۲۰۱۱) و سنتینل-۱ سری B‏ (دسامبر ۲۰۱۵ تا اوت ۲۰۱۷) بوده است؛ 
و نتیجه گرفتند که ۸۵ درصد از منطقه سیل‌زده با نرخ کمتر از ۵ میلی‌متر در سال فرومی‌نشیند. آن‌ها از طریق آزمون 
کای اسکوتر مستقل " بیان نمودند که فرونشست زمین از طریق تغییر پایه ارتفاع سیلاب" و گرادیان‌های شیب‌بر روی 
شدت سیل تأثیرگذار بوده است. حقیقی و معتق (۲۰۱۹) فرونشست دشت تهران و ورامین را مطالعه کردند. آن‌ها 
مطالعه جامعی را برای بررسی دقیق پدیده فرونشست به کمک روش ISAR‏ و پردازش YA‏ تصویر ماهواره انویست 
و ۱۰ تصویر ماهواره آلوس پالسار (ALOS-PALSAR)‏ و ۶۸ تصویر ماهواره تراسار-ایکس (TerASAR-X)‏ و ٦٤‏ 
تصویر ماهواره سنتینل-۱ برای یک بازه زمانی VE‏ ساله به کمک ۶0۰۰ اینترفرو گرام انجام داده‌اند. نتایج آن‌ها نشان داد 
که دشت تهران با ماکزیمم نرخ فرونشست YO‏ سانتیمتر در حال فرونشست است و فرودگاه بین‌المللی تهران سالانه 
۵ سانتی‌متر فرومی‌نشیند. ژئو" و همکاران (۲۰۱۹)» در مقاله‌ای به ارزیابی سری‌های زمانی الگوهای فرونشست زمین 
در دشت پکن شرقی» چین پرداخته‌اند. این محققین نخست. به روش تداخل سنجی پراکنش گر دائمی PSD‏ و با 
به‌کارگیری داده‌های سنجنده انویست و رادارست (Radarsat-2)‏ در بازه زمانی ۲۰۰۷ تا ۲۰۱۵ برای جابجایی‌های 


درازمدت این دشت استفاده گردیده است. نتایج این تحقیقی ماکزیمم نرخ فرونشست do‏ ۱/۱۸۶ میلی‌متر در سال تا 


1 Nikolaos Liosis 

2 Al Wagan 

3 Cristina Da Lio 

4 Miller 

5 The Chi-square Test of Independence 
6 Base Flood Elevations 

7 Junjie Zuo 
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سال ۲۰۱۵ به دست آوردند. [Jo‏ عمده الگوهای تکاملی سری زمانی فرونشست تغییرات سالانه تراز آب زیرزمینی 
بوده است. 

سیان ! و همکاران (۲۰۱۹)» فرونشست زمین در شهرهای ساحلی بانجول " (کامبیا) و اگوس " (نیجریه) در آفریقا 
را با استفاده از تداخل سنجی پراکنش گرهای دائمی (PSD)‏ و تصاویر سنتینل-۱؛ و تصاویر (TerraSAR-X)‏ تصاویر 
(COSMO-SkyMed)‏ با بهره‌گیری از پکیج خودکار (SNAP2ZSTAMPS)‏ را مورد مطالعه قرار داده‌اند. این روش 
خودکار برای مناطق با جابجایی زیاد زمین» اطلاعات مفید برای توسعه شهری, مدیریت خطر بلایا و برنامه‌ریزی 
سازگاری با آب‌وهوا را ارائه می‌دهد. فیاسچی * و همکاران (۲۰۱۹)» جابجایی زمین را با استفاده از روش پراکنش 
گرهای دائمی (PST)‏ در مناطق معتدل اقیانوسی با به‌کارگیری داده‌های راداری سنتینل-۱ در جمهوری ایرلند مطالعه 
کرده‌اند. بازه زمانی این مطالعه از سال ۲۰۱۵ تا ۲۰۱۸ بوده است. بیشترین مقدار جابجایی زمین حاصل از فرونشست 
را ۱۷- میلی‌متر در سال در مناطق شهری و محیطهای طبیعی رخ‌داده است. ال کمالی" و همکاران Y)‏ 
فرونشست و جابجایی زمین را در منطقه lazy‏ امارات متحده عربی به‌عنوان یک منطقه با آب‌وهوای SAH‏ و بیابانی 
بررسی کرده‌اند. سرعت فرونشست در این محدوده را ۶۰ میلی‌متر در سال به دست آورده‌اند که به علت افت ۱۲ 
متری سطح آب‌های زیرزمینی در این منطقه بوده است. در این تحقیق» از روش تداخل سنجی پراکنش گرهای پایدار 
(PS-InSAR)‏ بر مبنای پردازش خودکار ٩۸۳251۸۷۲۳5‏ برای برآورد فرونشست در دشت مرند استفاده شد. 

دشت مرند. یکی از دشت‌های مهم استان آذربایجان شرقی است که به دلیل بهره‌برداری بی‌رویه از آب‌های 
زیرزمینی جهت تأمین نیازهای آبی منطقه خصوصاً کشاورزی و همچنین کاهش نزولات جوی و خشک‌سالی‌های 
اخیر دچار فرونشست شده است. سطح آب‌های زیرزمینی در این دشت از سال ۱۳۹۰ تا AYAV‏ به‌طور متوسط تا Y‏ 
متر در سال برآورد شده است. که زمین‌های کشاورزی حاشیه کوشک سرای و یامچی بالاترین نرخ افت سطح 
ایستابی را داشته‌اند (داداشی و همکاران, ۱۳۹۹: (AY.‏ لذا با توجه به اهمیت این موضوع تحقیق حاضر با استفاده از 
۶ تصویر ماهواره ستتینل-۱ (Sentinel-1A)‏ برای به دست آوردن نرخ فرونشست زمین در محدوده مطالعاتی از 


سال ۲۰۱۹ تا ۲۰۲۰ انجام‌گرفته است که نتایج آن می‌تواند در مدیریت آبخوان این دشت مؤثر واقع شود. 


] Cian 

2 Banjul 

3 Lagos 

4 Simone Fiaschi 

5 Muhagir El Kamali 


سال یازدهم GL‏ پدیده فرونشست دشت‌ها بر مبنای الگوریتم خودکار .... ۳۷ 


۲- مواد و روش 
۱-۲- موقعیت منطقه موردمطالعه 


دشت مرند به‌عنوان بخشی از حوضه آبریز رودخانه ارس در شمال غرب ایران واقع‌شده است و ازنظر زمین 
ساختی به‌صورت یک گرابن - در حدفاصل رشته‌های کوهستانی قره Flo‏ در شمال و میشوداغ در جنوب - فروافتاده 
است. گسل‌های متعددی به‌ویژه گسل شمالی میشو در شکل‌گیری این دشت نقش داشته‌اند. این دشت با روند 
شرقی-غربی به‌صورت یک چاله ساختمانی نئوژن-کواترنری در شمال غرب ایران و حدفاصل دو رشته‌کوه عمده 
منطقه یعنی قره داغ و میشو داغ واقع‌شده است. این چاله انباشته از نهشته‌های آبرفتی دوران کواترنری است که در 
قالب مخروط افکنه هایی در امتداد جبهۀ کوهستانی مشرف به دشت ردیف شدهاند (مختاری AMAA ¢ SOS‏ 
قدیمی‌ترین واحد زمین‌شناسی منطقه» متعلق به پرکامبرین است که در جنوب دشت در مجاورت توده‌های نفوذی 
گرانیت شیل و گرانیت قوشچی از گسترش محدودی برخوردار است. دو رودخانه اصلی این زیر حوضه؛ زیلبر چای 
و زنوز چای است که پس از اتصال به رودخانه قطور چای و سپس ارس» به حوضه آبریز دریای خزر می‌ریزند. 
سلسله جبال جنوبی شامل ارتفاعات میشو و مورو هستند که حد بین حوضه‌های آبریز دشت مرند یکی از دشت‌های 
با توان UL‏ درزمینه تولید محصولات باغی است. این دشت بهره‌مندی مناسبی از منابع آبی داشته و این امس به رونق 
کشاورزی آن کمک نموده اسست. اما با توجه به مشکلاتی نظیر کم‌آبی و نیاز به افزایش تولید با توجه به افزایش 
جمعیت. مکانیزه شدن سیستم آبیاری» کاشت و برداشت محصولات در دشت مرند. در سال‌های اخیر دچار 


فرو 


شده است. برداشت‌های بی‌رویه از منابع آب زیرزمینی طی سنوات گذشته موجب کاهش ذخایر آبی و 


ایجاد مشکلات جدی ازجمله کاهش آبدهی locke‏ خشک شدن قنوات و چشمه‌هاء شور شدن آب‌های زیرزمینی در 


بخش‌های زیادی از CSS‏ مرند شده است. برداشت بیش‌ازحد آب و افت سطح آب زیرزمینی موجب ایجاد 3 dle»‏ 


و يا فرو زمین در این دشت شده است. فرو 


رخداده در دشت مرند باعث شده در منطقه بر 
زیرساخت‌های موجود تأثیر گذاشته و باعث شکست پایه بتنی چاه‌های کشاورزی و خم شدن تیرهای برق گردیده 


است. 
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شکل ۱- موقعیت جغرافیایی محدوده موردمطالعه. دشت مرند 


-Y-Y‏ روش انجام پژوهش 

در تکنیک پراکنش گرهای دائمی ou! (PSIBSAR)‏ اصلی شناسایی نقاط همبسته با رفتاری ثابت و پایدار در طول 
زمان با استفاده از سیگنال استخراج شده از داده‌های راداری است. درواقع بر روی هر پراکنش گر دائمی اطلاعاتی از 
قبیل نرخ سرعت جابجایی و خطای باقیمانده مدل ارتفاعی رقومی زمین برآورد می‌شسود روند کلی تکنیک پراکنش 
گرهای دائمی در شکل ۳ ارائه شده است. در مطالعه حاضر به‌منظور بررسی رفتار فرونشست در منطقه مطالعاتی» از 
پکیج StaMPS-PSI‏ استفاده شده است (هوپر» ۲۰۱۸). تغیبرات و افت سطح آب در طول YA‏ سال. از سال ۱۳۸۰ تا 
۸ با استفاده از داده‌های چاه‌های پیزومتری» از شرکت آب منطقه‌ای استان آذر بایجان شرقی دریافت و در سطح 
دشت بررسی و سطح آب در ۲۳ گمانه که داده‌های کامل و مشترک داشته‌اند در یک فایل اکسل تهیه گردید. سپس با 
استفاده از نرم‌افزار ArcGIS‏ روش درونیابی (DW)‏ انجام شد. از روش پراکنده ساز دائمی (PST)‏ بر مبنای پردازش 
خودکار SNAP2STAMPS‏ به‌منظور پایش فرونشست منطقه استفاده‌شده است. جهت پردازش از تصاویر ماهواره 
ستتینل-۱ مربوط به بازه زمانی ۵ ساله از تاریخ ۲۳ ژانویه ۲۰۱۹ تا ۲۲ دسامبر ۲۰۲۰ استفاده‌شده است. تعداد این 
تصاویر ۱۳۳ تصویر c pae‏ از نوع SLC‏ با قطبش (VV)‏ و حالت IW)‏ و پایین گذر (مسیر (VA‏ بوده است. بیس 
لاین زمانی هر تصویر. ۱۲ روز بوده است. به حداقل رساندن این بیس لاین به همدوسی داده‌ها می‌انجامد. شبکه 
گراف مربوط به تصاویر ستتینل-۱ هر پنج پارت تحقیق ترسیم شده است. گراف مورداستفاده در این پردازش به 
روش ستاره‌ای است (شکل Y‏ 


سال یازدهم پایش پدیده فرونشست دشت‌ها بر مبنای الگوریتم خودکار .... YA‏ 


شکل ۲-گراف ستاره‌ای اینترفروگرام های ایجادشده هرسال برای روش تداخل سنجی راداری (PSI)‏ 


در شکل ۲ نگاره‌های بیس لاين را برای هر پنج دوره موردبررسی در تحقیق را نشان می‌دهد. تصویر پایه در 
مرکز و به رنگ قرمز که بیشینه مقدار همبستگی زمانی» مکانی و نرخ داپلر را با سایر تصاویر دارد و موقعیت نسبی 
تصاویر پیرو را به رنگ آبی OLS‏ داده است. هر خط مخفف یک جفت تصویر است که برای محاسبه اینترفروگرام 
مورداستفاده قرار می‌گیرد. پردازش‌های لازم در نرم‌افزار اسنپ (SNAP)‏ و استمپس (STAMPS)‏ انجام گرفته است و 
ایتترفروگرام‌های لازم تهیه‌شده و درنهایت با استفاده از روش پراکنده ساز دائمی» میزان فرونشست به‌دست‌آمده است. 
مشخصات تصاویر انتخابی مورداستفاده در تحقیق در جدول Y‏ ارائه شده است: 


حدول ۱- مشخصات dS‏ تصاویر انتخایی مورد استفاده در تحفیق 


N-images | Fram Orbit Pass | Lastimage | Firstimage | Satellite sensors 
M 1٤ Descending v4 YV YYY ۱۲ ۱ YY Sentinel-1A 
YA £M Descending v4 YAN YrW/ASMMN Sentinel-1A 
Y. £M Descending v4 ۱ ۵ 11| Sentinel-1A 
۳۰ £M Descending v4 ۱ ۱/۰ AVT Sentinel-1A 
YA e Descending v4 YeY Y Yo Telefe) Sentinel-1A 


پردازش داده‌ها به روش پراکنده ساز دائمی در این تحقیق به دو بخش زیر تقسیم‌شده است که عبارت‌اند از: 
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پردازش خود کار با استفاده از یک تصویر پایه با استفاده از اسکریپت‌های بسته متن باز پایتون نرم افزار 


| ۰ ۱ 
لس 


پکیج منبع باز 5 توسط ژوزه مانوئل دلگادو بلاسکو" و میشائیل فیوملس "در سال ۲۰۱۸ 
dail!‏ است (دلگادو بلاسکو و همکاران ۲۰۱۹). در این تحقیق آخرین ورژن این پکیج که توسط دلگادو بلاسکو 
در سال ۲۰۱۹ ارتقاء یافت استفاده‌شده است (دلگادو بلاسکو و همکاران ۲۰۱۹). Ub py‏ پردازش خودکار 
5 طی چهار مرحله انجام می‌پذیرد. در این نوع پردازش» اسکرییت‌ها به‌طور سامانمند با استفاده از 
ابزار گراف» مجموعه‌ای از ایتترفروگرام را به حالت دسته‌ای * ایجاد می‌کنند. قبل از آن تصاویر پایه هر دوره انتخاب 
می‌شوند. انتخاب تصویر Master‏ به صورتی انجام می‌شود که در آن تأثیر عدم همبستگی‌های مکانی - زمانی و 
اختلاف داپلر مرکزی به حداقل برسد. آماده‌سازی و مرتب کردن تصاویر ستتیل-۱ بر اساس تاریخ LH‏ تصاویر انجام 
خواهند شد. بعدازاینکه تصویر Master‏ برای هر دوره انتخاب bo S‏ تمامی تصاویر Slave‏ بر اساس مختصات 
محدوده موردمطالعه برش داده‌شده و اطلاعات مداری آن‌ها به‌طور خودکار توسط SNAP‏ بارگیری و به‌روزرسانی 
من‌شوند S uem‏ 


TOPSAR-Split Apply-Orbit-File 


شکل ۳-گراف مربوط به مرحله جداسازی تصاویر و دریافت اطلاعات مداری آن‌ها در بسته متن باز پایتون نرم 


افزار اسنپ 


ثبت هندسی تصاویر لازمه تشکیل اینترفروگرام در مرحله بعدی است که تضمین می‌کند هر تارگت زمینی مربوط 
به یک پیکسل مشابه با مختصات آزیموت و رنج یکسان در هر دو تصویر پایه و پیرو است. مرحله بعدی ثبت 
هندسی داده‌ها و تولید اینترفروگرام ها است که با استفاده از تنوع طبقی آن‌ها صورت می‌گیرد (دلگادو بلاسکو و 
همکاران, ۲۰۱۹). این مرحله ay‏ لحاظ محاسباتی» دشوارترین مرحله پردازش در 831۸۳257۸۷1۳5 است. زیرا بت 
هندسی داده‌های تاپسار را انجام می‌دهد و هم‌زمان با تولید اینترفروگرام هاء مقدار سهم فایز نویز و فاز زمین مسطح 


نیز حذف‌شده و بخش‌های جداشده تصاویر نیز به همدیگر متصل می‌شوند (شکل 4). 


1 SNAP2STAMPS 
2 Delgado Blasco 

3 Michael Foumelis 
4 Batch mode 
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Master Read 
TopoPhase 
Removal 


Slaves Read 


شکل ۶- گراف مربوط مرحله ثبت هندسی و تولید اینترفروگرام ها در بسته متن‌باز پایتون نرم‌افزار اسنپ 


داده‌های تکمیلی موردنیاز استمیس ازجمله؛ Ub‏ ارتفاعی و مختصات عرض و طول جغرافیایی داده‌ها نیز تولید 


می‌شوند. درنهایت خروجی‌های استمپیس TIN‏ می‌شوند» که داده‌های پردازش‌شده مراحل قبلی را به فرمت 


رستری باینری و سازگار با استمپس تبدیل می‌کند (شکل ۵). 


شکل ۵-گراف مربوط مرحله خروجی گرفتن در بسته متن‌باز پایتون نرم‌افزار اسنپ 


StaMPS 
Export 


یکی از کاربردی‌ترین شاخص‌های شناسایی پیکسل‌های پراکنده ساز دائمی» شاخص پراکندگی دامنه هست که 
اش اویش eod‏ ارگ و اران وال uli Lal sat al E re‏ زا که palette‏ 
oe a oe TV] diu Giles diu yo a‏ کل اتکی راکد گی alb AA tulo‏ بیان BP sgh‏ 
و همکاران ۲۰۰۱). 


D, = % P] 


mA 


که در آن 504 م79 به ترتیب برابر با انحراف معیار و میانگین دامنه باز پراکنش هر پیکسل است. در این مرحله 
پراکنده سازهای کمتر ولی باقابلیت اطمینان بالا جهت برآورد اتمسفر انتخاب می‌شوند. مقدار بهینه این شاحص ۰/۶ 


پيشنهادشده aad‏ هر این تحقیق مقدار شاخصی براکندگی دامنه EY)‏ ) انتخاب‌شده استه تا تعداد براکنده سال‌های 


] Ferretti 


yy‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


دائمی) اولیه را در آغاز پردازش افزايش دهد این دامنه به صسورتی تعیین می‌شسود که احتمال انتخاب پراکنش 
کننده‌های واقعی را افزايش داده و احتمال انتخاب پیکسل‌های با فاز تصادفی کاهش Ll‏ مرحله بعد از آماده‌سازی 
داده‌های استمیس, مربوط به واردکردن داده‌های آماده‌شده در اسنپ به استمپس است. پس‌ازآن» زنجیره پردازش 
پراکنده ساز استمپس (StaMPS-PSI)‏ ایجاد می‌شود (هوپر" و همکاران» ۲۰۱۸). 

پردازش روش پراکنده ساز دائمی با استفاده از استمپس (StaMPS)‏ 

مجموعه‌ای از اینترفرو گرام ها با استفاده از روش‌های اجرایی متلب در استمپس پردازش‌شسده اسست. که با 
تصحیحات تروپوسفریک مبتنی بر فاز حطی با استفاده از پکیج ترین TRAIN)‏ که تأخیر تروپوسفری را نسبت به 
تویوگرافی جبران می‌کند. تکمیل گردیده است (بکرت" و همکاران» ۲۰۱۵). 

مراحل پردازش تصاویر راداری سنتینل-۱ در c‏ استمپس Gb StaMPS)‏ شکل A‏ شرح داده شده است: 


A پردازش سری زمامی‎ N, S لاح فار‎ P 
داده ها در‎ ‘a 
s A فاز‎ gm = 
a —7 p V 
——— 71 e بر آورد خطای ;444 دید‎ 2 
es تخسین نوبز فا‎ » — 
م‎ at) های‎ PS Usd ™ " Master »" 


براورد daii‏ جایجاب 


She PN <á " T 
(StaMPS) مراحل پردازش استمپس‎ -٦ شکل‎ 
نتایج و بحث‎ -Y 
iaa به‌منظور بررسی نوسانات سطح آب زیرزمینی در این محدوده. با استفاده از اطلاعات ثبت‌شده سطح آب‎ 
نقشه هم افت‎ (AYAY در ۳۰ چاه مشاهداتی که دارای داده‌های کامل بوده است. در بازه زمانی ۷ ساله (۱۳۸۰ تا‎ 
درون‌پابی به روش ۷ ترسیم‌شده‎ iub و از‎ ARC-GIS نرم‌افزار‎ daa سطح آب زیرزمینی محدوده مطالعاتی در‎ 


طی 0592 V‏ ساله بوده sul‏ که ازنظر Nm‏ مکانی. مناطق میانی دارای بالاترین میزان افت بوده‌اند. جاه‌های واقع 


] Hooper 
2 Bekert 


سال يازدهم پایش پدیده فرونشست دشت‌ها بر مبنای الگوریتم ET‏ 


در محدوده‌های شمال غرب و جنوب غرب کشک سرای و یامچی دارای بالاترین افت سالانه آب زیرزمینی را داشته 


است. 


| آب سالانه امتر نقاط شهری‎ cil 


$ d py iGo 
<< جاده اصلی‎ TTET 


شکل ۷- نقشه افت منابع آب زیرزمینی منطقه موردمطالعه 


همچنین برای بررسی تغییرات زمانی نوسانات سالانه سطح آب زیرزمینی از داده‌های ماهانه سطح آب زیرزمینی 
ole ۰‏ مشاهده‌ای در طی سال‌های ۱۳۸۰ تا ۱۳۹۷ استفاده‌شده است. شکل A‏ هیدر و گراف مربوط به تغییرات سری 
زمانی سطح آب در ۳۰ چاه مشاهده‌ای موجود در دشت مرند را نشان می‌دهد. بررسی هیدروگراف‌های چاه‌های 
مشاهده‌ای قرارگرفته در آبخوان مرند نشان‌دهنده حالت نوسانی تراز آب چاه‌ها است و تراز آب همواره روند نزولی 
نسبت به زمان نشان می‌دهد. روند نزولی هیدروگراف ها نشان‌دهنده افت مداوم سطح آب زیرزمینی است. کاهش 
تغذیه ناشی از خشک‌سالی‌های اخیر به همراه افزایش برداشت بی‌رویه به‌عنوان عامل اصلی افت شدید سطح آب 


زیرزمینی در دشت مرند شده است. 


شمارة چهارم 


جغرافیا و مخاطرات محیطی 


yé 


POTET 
LM 
N- 
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۱ 4 
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۳ 
T 
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شکل ۸- نمودارهای سری زمانی چاه‌های مشاهداتی دارای افت آب و فرونشست در محدوده مطالعاتی 


تحلیل نقشه کاربری و پوشش اراضی و ارتباط ان با فرونشست منطقه 


نوع کاربری اراضی از عوامل انسانی مهم در تشدید وقوع فرونشست است. درواقع» بهره‌برداری‌های انسان یکی 


از عوامل موثر در وقوع فرونشست محسوب می‌شود. به‌طوریکه مناطق دارای کاربری اراضی کشاورزی آبی به دلیل 


سال یازدهم پایش پدیده فرونشست دشت‌ها بر مبنای الگوریتم Yo ETT‏ 


بهره‌برداری از مناطق آب زیرزمینی» پتانسیل بالایی جهت وقوع فرونشست دارد. همچنین مناطق سکونتگاهی نیز به 
دلیل ساخت‌وسازهاء مستعد وقوع فرونشست هستند. در این تحقیق, نقشه کاربری اراضی محدوده موردمطالعه از 
s a‏ اراس کل ag CLS oos aras oes‏ که ۲۷۲۱ از اک با Soto‏ کر gall pas‏ 
داده‌های راداری Sentinel-1‏ و Sentinel-2‏ در ۱۱ کلاس پوشش زمین تهیه‌شده اخحذ گردیده و محدوده مطالعاتی از 
تشه fs‏ پراش داد کد is‏ موی Jr‏ 258« سومان BSA‏ ساسا ریش کے lus outs‏ 
ARC-GIS‏ از طریق ابزار Reclassify‏ به ۵ طبقه نواحی سکونتگاهی: کشاورزی - باغی. مراتع و اراضی بایر 
طبقه‌بندی‌شده است (شکل 4( بر اساس نقشه مذکون مناطق جنوبی و جنوب غربی محدوده منطبق بر کوهستان 
میشو را کاربری مراتع و اراضی بایر و همچنین مناطق مرکزی محدوده را اراضی کشاورزی آبی و ديم دربرگرفته 
است. فرونشست رخداده در این منطقه در کلاس مربوط به اراضی زراعی نزدیک به مناطق سکونتگاهی اتفاق افتاده 
است. با توجه به این تغییرات» میزان افت آب‌های زیرزمینی و یا میزان افزایش of‏ را مشخص سازد. مطالعات 
بیشتری در کل جهان و کاربری‌های مختلف تأثیرات بسیار زیادی بر روی سطح تراز GET‏ زیرزمینی دارند. 


به این دلیل که کاربری‌های مختلف باعث مصرف متفاوتی از آب‌های زیرزمینی می‌شود. 


شکل -s‏ نقشه کاربری اراضی محدوده مطالعاتی 


۳۹ جغرافیا و مخاطرات محیطی شمارة چهارم 


ارزیابی ON x‏ فرونشست L‏ استفاده از روش تداخل E ox‏ گرهای دائمی (PSInSAR)‏ 


بررسی وضعیت فرونشست در محدوده موردمطالعه در بازه زمانی VY‏ ژانویه سال ۲۰۱ تا ۲۳ دسامبر سال ۲۰۱ 


نشان می‌دهد که بیش‌ترین میزان 2339 ت در قسمت میانی دشت مرند. در قاعده مخروط افکنه های قدیمی و 


قسمت دشتی این محدوده. در محدوده روستاهای اربطان, فارفان AS yo‏ یالقوزآغاج. گله بان و درویش محمد بوده 
است. حداکثر فرونشست در سال Y* YA‏ با فاصله زمانی ۲ روزه ۱۳۷/۸- میلی‌متر در سال برآورد شد. کمترین میزان 


جابجایی با مقدار ۱۹/۱ میلی‌متر بالاآمدگی در منطقه موردمطالعه رخ‌داده است (شکل ۱۰). 


19.1 


- 137.8 


45.3 45.4 455 45.6 45.7 45.8 45.9 46 


شکل ۰- نقشه سرعت حابحایی متوسط محدوده مطالعاتی سال ۲۰۱۳۲۱ 


نقشه g‏ 23 = سال ۷ حاکی از افزايش نرخ فرونشست در محدوده مطالعاتی است (شکل ۱ حداکثر 
جابجایی زمین در سال ۲۰۱۷ در بازه زمانی ۲ روزه. ۲ - میلی‌متر به‌صورت فرونشست و ۲۵ میلی‌متر به شکل 
بالاآمدگی بوده است. تغییرات جابجایی زمین در این سال نسبت به سال قبل بدین‌صورت بوده است که نرخ 


فرونشست ۱۸ میلی‌متر کاهش و بالاآمدگی منطقه " میلی‌متر افزایش پیداکرده است. 


سال یازدهم پایش پدیده فرونشست دشت‌ها بر مبنای الگوریتم خودکار .... YV‏ 


25 


mm/yr 


-120.2 


454 4*3 456 i57 45 8 159 16 161 at 


شکل ۱- نقشه سرعت حابحایی متوسط محدوده مطالعاتی سال ۲۰۱۷ 


تاریخ ۱۱ ژانویه ۲۰۱۸ تا ۲۵ دسامبر ۲۱۰۸ میزان فرونشست منطقه ۳۲ میلی‌متر افزایش‌یافته و به نرخ حداکثری 


۲ - میلی‌متر فرونشست و ۳۶/۲ میلی‌متر بالاآمدگی رسیده است AVY JS)‏ 


31.5 


mm/yr 


3445 


-150.2 


شکل ۲- نقشه سرعت حابحایی متوسط محدوده مطالعاتی سال ۲۰۱۸ 


بررسی وضعیت فرونشست در محدوده موردمطالعه از روز V‏ ژانویه سال ۲۹ نا ۰ دسامبر همین سال نشان 


می‌دهد که نرخ جابجایی به‌صورت /۱۲۳- میلی‌متر فرونشست و ۲۲/۵ میلی‌متر بالاآمدگی به‌دست‌آمده است؛ 


وضعیت d A53‏ زمین در این محدوده روند افزایشی NY‏ میلی‌متر بالاآمدگی و ۲۹ میلی‌متر کاهش فرونشست 
زمین نسبت به سال قبل را نشان می‌دهد (شکل ۱۳). 


۳۸ جغرافیا و مخاطرات محیطی شمارة چهارم 


22.6 


mm/yr 


-123.6 


شکل ۱۳-نقشه سرعت جابجایی متوسط محدوده مطالعاتی سال ۲۰۱۹ 


پررسی وضعیت فرونشست در محدوده موردمطالعه در بازه زمانی روز Y‏ ژانویه سال ۲۰۲۰ تا دسامبر سال ۲۰۲۰ 
نشان می‌دهد که این مقدار نسبت به سال قبل حدود ٩‏ میلی‌متر کاهش فرونشست و ۱ میلی‌متر بالاآمدگی افزايش 
پیداکرده است. در سال ۲۰۲۰ مقدار فرونشست و بالاآمدگی به ترتیب ۱۳۱/۷- میلی‌متر و ۲۳/۱ میلی‌متر به‌دست‌آمده 
است (شکل OE‏ بیش‌ترین میزان جابجایی رخ‌داده در محدوده مطالعاتی» در سال ۲۰۱۸ اتفاق افتاده است درمجموع 


مخروط افکنه‌ها Le par‏ قاعده مخروط افکنه‌ها و دشت انتهایی دارای بیش‌ترین مقادیر فرونشست هستند. 


mm/yr 


-131.7 


شکل ۶- نفشه سرعت حابحایی متو سط محدوده مطالعاتی سال ۲۰۲۰ 


سال یازدهم GL‏ پدیده فرونشست دشت‌ها بر مبنای الگوریتم خودکار .... ۳۹ 


۶- جمع‌بندی 

روش سری زمانی تداخل سنجی با پراکنش گرهای دائمی (PSD)‏ کمک شایانی در جهت پایش فرونشست و 
جابجایی‌های سطح زمین کرده است. ازجمله محدودیت‌های این روش وجود پیکسل‌های سفید زنگ در سطح نقشه 
جابجایی بود. این موضوع به دلیل ناهمبستگی زمانی که براثر تغیبرات سریع پوشش منطقه و درنتیجه پایین بودن 
میزان همدوسی ایجاد می‌شود. همچنین اختلاف‌فاز ناشی از باز کردن (phase unwrapping) jb‏ که با توجه به 
توپوگرافی منطقه موردمطالعه قابل‌توجه نبوده است. در تحقیق حاضر برای نخستین بار به پایش نرخ فرونششست 
زمین در دشت مرند استان آذربایجان شرقی» از فن تداخل سنجی راداری PSI‏ با استفاده از داده‌های سنجنده سنتینل- 
| و پکیج خودکار SNAP2STAMPS‏ استفاده‌شده است. همچنین در این تحقیق, پتانسیل بسته متن‌باز استمپس و 
نرم‌افزار اسنپ جهت پردازش تداخل سنجی راداری OLS‏ داده‌شده است که از نتایج مهم استخراج‌شده از این 
الگوریتم پردازشی, ایجاد نقشه‌های متوسط جابجایی برای دشت مرند با مقادیر NY AYN‏ ۱۲/۳-۱۵/۲ - و 
- سانتی‌متر در سال در راستای دید ماهواره در طول بازه زمانی موردمطالعه (از سال ژانویه Y VV‏ تا دسامبر 
۰ بوده است. فرونشست‌های رخ‌داده در منطقه بر زیرساخت‌های موجود تأثیر گذاشته و باعث شکست A‏ بتنی 
چاه‌های کشاورزی و خم شدن تیرهای برق گردیده است. به لحاظ پراکنش مکانی, مناطق فرونشست یافته عمدتاً در 
مناطق میانی دشت مشاهده‌شده است؛ بنابراین مناطق میانی دشت بیشترین نرخ فرونشست نسبت به سایر نواحی 
دشت را دارا هستند. منطقه فرونشست در جنوب شهر یامچی و شمال و شمال غرب کشک سرای در حال اتفاق 
افتادن است و روستاهای زیادی ازجمله اربطان» فارفارء مرکید ایلات یالقوزاغاج» درویش محمد در آن واقع‌شده در 
این کاسه فرونشست قرار دارد که منطبق بر کاربری اراضی کشاورزی و باغی است که به سمت مناطق شهری در 
حال پیشروی است. بررسی هیدروگراف‌های چاه‌های مشاهده‌ای قرارگرفته در آبخوان مرند نشان‌دهنده حالت نوسانی 
تراز آب چاه‌ها است و تراز آب همواره روند نزولی نسبت به زمان نشان می‌دهد. که نشان‌دهنده افت مداوم سطح 
آب زیرزمینی است. فرونشست‌های رخ‌داده در دشت مرند بیشتر منطبق بر زمین‌های کشاورزی بوده و سرعت بالایی 


دارد. 
کتاینامه 
احمدی, سلمان؛ سودمند افشار» رضا؛ ۱۳۹۹. پایش فرونشست دشت‌های قروه و چهاردولی استان‌های همدان و کردستان 


به دلیل برداشت بی‌رویه از منابع آب‌های زیرزمینی با استفاده از فن پراکنش گرهای دائمی. محیط‌زیست و مهندسی 
آب. دوره A‏ شماره Y‏ صص ۲۱۹-۲۳۳ 
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ژئوتکنیک. مشهد. صص Met‏ 
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سنجی راداری. پایان‌نامه کارشناسی y‏ دانشگاه صنعتی خواجه‌نصیرالدین طوسی. دانشکده مهندسی نقشه‌برداری 
(زئودزی و ژئوماتیک). صص AVe‏ 
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شماره Y‏ صص ۸۵-۹۱۰ 

شرافت. متین؛ انصاری, عبدالحمید؛ مجتهد coals‏ سید حسین؛ قربانی احمد؛ MAS‏ پایش فرونشست دشت ابرکوه یزد 
با استفاده از فن تداخل سنجی راداری مبنی بر پراکنش گرهای LL‏ نحشک بوم. دوره A‏ شماره ۲. صص AYIYI‏ 

شمشکی امیر؛ انتظام سلطانی. ایمان؛ بلورچیء محمدجواد؛ ۱۳۸۶ فرونشست زمین در دشت تهران و عوامل موثر در 
شک لگیر ی آن.. سازمان زمین‌شناسی و اکتشافات معذنی کشور. ص ۲۱ 

علی بخشی» حمید؛ ۱۳۹۵. بررسی میزان فرونشست دشت ورامین با استفاده از ابزارهای ژئودتیکی داده‌های ترازیابی 
دقیق NSAR‏ و تصاویر راداری GPS‏ پایان‌نامه کارشناسی ارشد. دانشگاه آزاد اسلامی. واحد شاهرود. دانشکده 
فنی و مهندسی گروه عمران. گروه ژئودزی. ص ۱۳۱. 
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